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CHAPTER 1 


INTRODUCTION 

Prediction of airplane response to atmospheric turbulence is 
commonly done with the assumption that the turbulence velocity varies 
along the flight path but does not vary along the span. In some previous 
studies where spanwise turbulence has been considered, power spectral 
methods were generally used (Coupry 1972, Fuller 1968). However, this 
investigation of the effects of spanwise turbulence is carried out in a 
different manner. Based on flight data, the aerodynamice forces and 
moments on the wing of NASA's B-57 aircraft are calculated using a modi- 
fied vortex strip theory. This model allows the input airspeed and 
angle of attack to vary along the wing. Flight data provide three 
spatially separated records of the two input parameters. As a first 
order approximation, linear distributions in angle of attack and air- 
speed are assumed between the three recording locations; at the nose 
and at each wing tip (see Figure 1.1). For details on the B-57 instru- 
mentation see Camp et al . (1984), Painter and Camp (1983), and Campbell 
et al. (1983). The modified vortex theory, hereafter referred to as 
the strip theory, has the ability to calculate lift and induced drag as 
a function of spanwise position. Furthermore, as in situations where 
these distributions are asymmetric, see Figure 1.2, an integration of 
these two forces over the wing span yields roll and yaw moments acting 
on the wing solely due to the asymmetric distribution of lift and drag. 
The results from strip theory calculations are next incorporated into a 
six-degrees-of-freedom (6D0F) program to investigate the effect of 

spanwise turbulence on aircraft response. 
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Figure 1.1. Spanwise variation of relative wind and angle of attack. 
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Figure 1.2. Asymmetric distributions of lift and drag. 


Following, the procedure of introducing the strip theory calculations 
into the 6D0F simulation is explained in more detail. From 40 seconds 
of flight data lift, induced drag, and roll and yaw moment time histories 
are calculated based upon 0.25-second averages of the input parameters. 

The time histories of roll and yaw moments represent moments acting on 
the aircraft wing that would be unaccounted for if no spanwise variation 
in turbulence were considered. These moments are referred to as addi- 
tional, and the roll moment expressed as: 

s 

r 


l add 


^(y)ydy 


(l.l) 


-S 


Similarly, the yaw moment is given by: 
s 


^ADD 


d i (y)ydy 


-S 


J 


( 1 . 2 ) 
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where s is the half-span length and y is the lateral axis in the body 
axis coordinate system. The &(y) and d^y) indicate sectional lift and 
induced drag as a function of y. 

Precalculated time histories of these two moments are then added 
into the force and moment equations in the 6D0F program (see Chapter 
IV). It is assumed that wing pitching moment is unaffected by the span- 
wise distribution of aerodynamic forces.* Furthermore, additional lift 
and induced drag are added as spanwise turbulence effects. The additional 
lift is given by: 

L ADD = L SPANWISE " l uniform 
and induced drag is expressed: 

D- = D i - D 1 (1 .4) 

1 ADD 1 SPANWISE 1 UNIFORM 

The subscript SPANWISE indicates total quantity of the aerodynamic force 
as calculated using strip theory with recorded variation in angle of 
attack and airspeed as input. The uniform terms come from a second 
calculation using strip theory, but this time the three separate measure- 
ments of the input parameters are averaged at each time step. These 
averages are then input to the strip theory model giving total lift and 
induced drag in a spanwise uniform wind field. Equations 1.3 and 1.4 
thus are expressing the effect spanwise turbulence has on lift and induced 
drag . 


*This is justified considering the negligible moment arm between 
aircraft center of gravity and the wing aerodynamic center. Changes in 
total lift due to spanwise turbulence are also small. 
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In Chapter II of this report the derivation of the strip theory is 
carried out. Chapter III contains the results of applying the strip 
theory to data from two test flights. The procedure for the 6D0F simu- 
lation of the test flights is described in Chapter IV. Results presented 
in Chapter IV demonstrate how spanwise turbulence effects influence the 
flight simulation of the B-57 aircraft. The last chapter contains the 
conclusions of this study. 
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CHAPTER II 


DERIVATION OF STRIP THEORY 

Vortex theory is the method most commonly used to calculate the 
lift and drag distributions of a finite wing. The strip theory developed 
in this chapter is a modified vortex theory which allows for random 
spanwise variations in angle of attack and airspeed. As the name implies, 
the airfoil is divided into a finite number of strips. At each strip, 
lift and drag can be evaluated. By integrating these forces across the 
wing span, the desired moments are obtained. 

1. DERIVATION OF TOTAL QUANTITIES 

The lift distribution across a wing can be represented by a distri- 
bution in the strength of the bound circulation r(y) as shown in Figure 
2.1. At any position in the spanwise direction, the lift per unit span 
can be expressed by Zhukovsky's theorem as (Kuethe and Chow 1976): 

*(y) = pV(y)r(y) (2.1) 

where p is the air density and V (y ) is the speed of the wing along the 
direction of the x-axis in the wind coordinates. The frame of reference 
in wind coordinates is moving with the mean wind. The x-direction is 
chosen to be parallel to the direction of flight. Contrary to tradi- 
tional analysis, V (y ) is considered a variable in the body axis y- 
direction. Similarly the induced (vortex) drag per unit span at any 
location again by Zhukovsky's theorem becomes: 

d -j (y) - pw(y)r(y) (2.2) 
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Figure 2.1. The relation between spanwise load variation and trailing 
vortex strength. 
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where w(y) is the downwash velocity. The total lift is found by inte- 
grating over the span: 
s 

L = f pV(y)r(y)dy (2.3) 

-s 

In the same manner, total induced drag is: 
s 

D i = pw(y)r(y)dy (2.4) 

« 

-s 

As a first approximation, the circulation is often assumed to be 
elliptical in nature, which leads to closed form solutions of the inte- 
grals. In this study the integrals are solved numerically since the 
circulation can be non-elliptical and asymmetric. 

The lift per unit span as in Equation 2.1 is also written: 

*(y) = \ pV 2 (y)c(y)C L (y) (2.5) 

where c(y) is the chord length and C L (y) is the sectional lift coeffi- 
cient. Combining Equations 2.1 and 2.5 and solving for r(y), the 
following expression is obtained: 

r(y) = \ V(y)c(y)C L (y) (2.6) 

The continuation of the derivation is simplified by designating the 
variable y in terms of the angle e as illustrated in Figure 2.2. 

Equation 2.6 is now written: 

r(e) = 4sV ( 0 ) p sha p e ( 0 ) ( 2 - 7 ) 

where p shape is a dimensionless shape parameter which includes the 
variation in both lift coefficient C L (y) and the chord length c(y). The 

circulation is made proportional to four times the half-span length 
rather than the root chord length for convenience. 

i 8 
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The shape parameter is represented as a Fourier sine series: 

p shape (e) * j, sin ne (2 ’ 3) 

The sine series satisfies the end conditions of the curve reducing to 
zero at the tips where ne is zero or n (n is an integer). A n is the 
amplitude of each sine curve making up the shape parameter. In Figure 
2.3 a sample asymmetrical shape parameter is illustrated. 

The total quantities of lift, induced drag, and roll and yaw 
moments given by Equations 2.3, 2.4, 1.1, and 1.2, respectively, are now 
written by substitution of Equations 2.7 and 2.8 as: 


L = l A 

n=l r 


4s 2 pV 2 (e)(sin ne sin e)de 


(2.9) 


.3 


9 



p 

shape 



D i 


oo 



00 

’ 2 

4s pV(e) (sin ne sin e)w(e)de 

i 

0 


o O 

£ Ann = I A n I 4s pv (e)(sin ne sin e cos e)de 


n=l 


if 

00 f 3 

^ADD = £ A n 4s pV(e)(sin ne sin e cos e)w(e)de 
n=1 L 


The change in variable going from y to e is presented below: 

0 7T 


S 0 

dy = 


dei 


de = 


(-s)sin e de = 


s sin e de 


-s 


(2.10) 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 
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The obvious tasks then becomes to find solutions for the A coeffi- 

n 

cients as well as to evaluate the integrals in Equations 2.9 through 
2.12. As already mentioned, the approach taken is to divide the wing 
into a finite number M sections (Figure 2.4). The integrals are then 
evaluated numerically using Simpson's rule. As will be explained in the 
following section, the number of A n coefficients that can be solved for 
is the same as the number of strips, M. Thus, the Fourier series repre- 


sentation of P shape 


is approximated to the desired accuracy by choosing 


a sufficiently large number M. 



2. DETERMINATION OF THE A n COEFFICIENTS AND DOWNWASH 

The lift coefficient at a section y along the span in terms of the 
lift curve slope is given as: 

C L (y> = aJCc,(y) - % (y!) - sCy)] (2.14) 

where a is the two-dimensional lift curve slope; a(y) is the angle of 

oo 

attack; a Q (y) is the section zero lift angle; and e(y) is the downwash 
angle. Since wings commonly are designed with washout or a variable 
incidence angle, a Q (y) is written as a variable of y. The angle a Q (y) 
is negative and adds to the measured angle of attack a(y). The downwash 
angle e(y) is positive and reduces the effective angle of attack a eff (y). 
The geometry of Equation 2.14 can be visualized by use of Figure 2.5. 

Solving Equation 2.6 in terms of C L (y) and equating it with Equa- 
tion 2.14 leads to the following result: 

V(y)[a(y) - a o (y)] = + w(y) (2.15) 



Figure 2.5. Effective angle of attack. 
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The downwash angle e(y) is given by: 


e(y) = tan -1 (2.16) 

Using a small angle approximation downwash becomes: 

w(y) = V(y)e(y) (2.17) 

The downwash term in Equation 2.15 must be expressed as a function of 
circulation before solutions to the A n coefficients can be found. Using 
the Biot-Savart law, velocity induced at a point y due to a semi- 
infinite vortex line is given by: 



where r is the circulation strength and h is the perpendicular distance 
from the vortex line to the point y (Figure 2.6). Consider now the 
influence on downwash at a position y along the span of the wing (Figure 
2.1, page 7) due to a trailing vortex filament of strength $r shed at 
y-| • Using Equation 2.18 this is expressed: 
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(2.19) 


dr(y,) 

TyT-^l 


6w(y) ' 4„(y’ - y) 


Since all trailing vortices influence the downwash velocity at position 
y, the total downwash is found by integrating Equation 2.19 across the 
span: 

drCyJ 

1 f 


dy, dy l 


-s 


y] - y 


( 2 . 20 ) 


Changing the variable and integrating Equation 2.20 by parts the down- 
wash becomes. 


w(e) = - 


1 


it V ( 0 1 ) l A n sin ne-| sin de-j 


( 2 . 21 ) 


(cos e 1 - cos 0 ) 

Finally, by substituting Equation 2.21 into Equation 2.15 and expressing 
the circulation in terms of the Fourier series, the following equation 
is obtained: 


p(e)[a(e) - a (e)] * 


n=l A|1 


sin ne - 


v(0) J sin ne l sin 9 l 

:v ‘ 9 ) J (COS 9, - cos e) 2 d61 


( 2 . 22 ) 


where n(e) is: 

c(e)a m 


>( 0 ) = 


8s 


(2.23) 


At this point in the derivation the finite number (M) strips are intro- 
duced. Assuming that V(e), c(e), a(e), and a Q (e) all are known at each 
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spanwise position m, Equation 2.22 is rewritten as: 


u(e m )[a(e m ) - a (ej] = T A n sin 
' nr ' m o' nr n 

n= I 


y( 9 m) 7 V(0,) sin 0 , sin n0, 

Hi » ! L Hft 

q (cos 0^ - COS 0 m ) 2 \ 

Solution to the A n coefficients is then given by: 


(2.24) 


M -1 

■ I v(9 m )[a(e m ) - a 0 ( e m )^nm 


(2.25) 


where I nm is the matrix identifiable in Equation 2.24. Immediately, 
solutions of Equation 2.25 only exist when N equals M. Furthermore, the 
singularity inside the integral at 0-| = 0 m must be avoided. For most 
convenient numeric evaluation, I nni is written: 


'nm = s1n ne m " “Tire 


n si n ne„ 


1 ] (V(e 1 ) - V(0 m )) sin 0 ] sin n0 ] 




(cos e.| - cos 0 m )‘ 


(2.26) 


The downwash given by Equation 2.21 can now be calculated at each of the 


spanwise sections as follows: 


S . v ‘ e n,> 


• n I, % 


[I rm - s1n ne m ] 


(2.27) 


3. COMPUTATION OF TOTAL QUANTITIES 


Equations 2.11 and 2.12 can now be solved for and N^, 

respectively. The quantities inside the integrals are all known at M 
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finite positions along the span. The integrals are then estimated 
numerically and the summations are performed with a finite number N A n 
coefficients. Equations 2.9 and 2.10 are solved in a similar manner. 

The two solution sets, the first obtained using a spanwise distribution 
of airspeed and angle of attack and the second obtained using a uni- 
form distribution, are then input into Equations 1.3 and 1.4, and the 
additional lift and drag forces are calculated. 

In this study the wind is divided into 21 strips. Thus, the 
computations involved, in addition to estimating integrals numerically, 
include a 21 by 21 matrix inversion. A computer code written in Fortran 
has been developed to execute the large number of calculations. 



CHAPTER III 


STRIP THEORY CALCULATIONS 

The reason for developing a strip theory with angles of attack and 
airspeed as spanwise variables is to utilize the data that has become 
available from the NASA Gust Gradient Program. Using the theory developed 
in the previous chapter, calculations of spanwise turbulence effects can 
now be performed. In this chapter the computation procedures are 
described and results presented. 

1. DATA AVAILABLE 

The NASA B-57 research aircraft is equipped to measure a large 
number of quantities. Table 3.1 lists all the variables directly 
recorded or calculated in the format received. Of interest in this 
part of the study are the angles of attack and airspeeds. The B-57 
is capable of recording these values at three different positions. This 
is made possible by installing flow vanes and pitot tubes at each wing 
tip in addition to the nose (see Figure 1.1, page 2). 

The variables are recorded on analog tapes at a rate of 200 per 
second aboard the airplane. Then at NASA Langley Research Center, 

Hampton, Virginia, digitized magnetic data tapes in engineering units 
are produced containing 40 sample points per second. In this study a 
10-point averaging scheme is used resulting in four samples per second. 

Each test flight consists of a number of runs. An individual run 
represents a continuous period of the flight at which data was recorded. 
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The length of a run may vary between less than half a minute and several 
minutes. 

This first part of the study analyzes two runs originating from two 
separate flights. They have both been chosen because of high levels of 
turbulence. The first is Run 21 of Flight 6 recorded July 14, 1982, at 
Denver, Colorado, as part of the Joint Airport Weather Studies (JAWS) 
Project. The other is Run 18 of Flight 60 recorded February 1, 1984, at 
Boulder, Colorado, during the Orographic Effects Campaign. Additional 
information on the flights is listed in Appendix D. In these two runs 
the aircraft flew level flights approximately 1000 feet above the 
ground. The data tapes from Flight 60 include records of pilot control 
inputs but Flight 6 does not. For the response analysis in Chapter IV, 
this is of significance. 


2. CALCULATION DESCRIPTION 

The strip theory developed in Chapter II has been programmed in a 
computer code referred to as Program I. Inputs to this program are data 
files from the B-57 test flights containing time histories of angles of 
attack and airspeed at the three different spanwise positions. At 
each time step of the input file the strip theory calculations are 
performed. No aircraft dynamics enter into Program I; each time step is 
simply treated uniquely, neglecting unsteady effects. The output file 
is a time history of lift, induced drag, and roll and yaw moments equal 
in length to the input file. The rate at which calculations can be done 
on the Pixel Supermicro Motorola 68,000computer is substantially slower 
than required to keep up with inputs in real time. Each time step of 
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0.25 seconds will consume approximately one minute of computer time. 

This is largely due to the inversion of the sizable 21 by 21 matrix. 

The results of Program I are stored, and as described in Chapter IV, a 
dynamic response analysis based on these results is carried out. 

The calculations in this chapter are designed so that the results 
can be directly utilized in FWG's six-degrees-of-freedom performance 
analysis computer program. The input to the 6D0F program representing 
spanwise turbulence is the quantities defined by Equations 1.1 through 
1.4. Two iterations using the strip theory are necessary. In the first, 
the recorded linearized distribution of angle of attack and airspeed at 
each time step are input, see Figure 3.1. The first iteration yield time 

histones of L m> ^ add > l $PANWISE’ and D i spANWISE ’ Next ’ 1n a second 
iteration, the calculations are based on inputs uniform across the span 
representing the average of the three spanwise recordings (Figure 3.2) 
uniquely determined at each time step. The second calculation gives 
L-iiuTFnDM and D. but no moments as they vanish when the distribution 

'uniform 

in lift and drag are symmetric. The two iterations using the strip 
theory then give the time histories of the four quantities representing 
spanwise turbulence needed for the response analysis, 

3. ANALYSIS OF STRIP THEORY RESULTS 

In the analysis of strip theory results, two problems were uncovered. 
First, the value of lift calculations appeared to be exaggerated as a 
result of neglecting unsteady effects. Secondly, the angles of attack 
seemed to contain a constant offset. The methods used to solve these 
complications are presented preceding the final results in this section. 
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Left Tip 



Figure 3.1. Sample of a linear distribution of inputs to Program I. 



Spanwise airspeed distribution Uniform airspeed distribution 

Figure 3.2. Spatial and uniform spanwise distributions of airspeed. 
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At each step of the calculations in Program I a steady-state 
condition is assumed. It would then be expected that the instantaneous 
lift and drag values are overestimated in gusts. A short literature 
review led to the estimation of a gust alleviation factor K, accounting 
for unsteady effects. This factor K = 0.7 is applied to the results as 
described in Appendix A. 

A second complication attributed to the inclination angles of the 
probes on the B-57 introduces offsets in the recorded angles of attack. 
Because of unrealistic results from the calculations applied to Run 18 
of Flight 60 an investigation was initiated. It showed that the time- 
averaged angle of attack for Flight 60, Run 18, differs between the three 
recording probes. Specifically, the left probe indicates an angle of 
attack about two degrees higher than the center probe, whereas the right 
probe lies close to half of a degree below on a time-averaged basis. 

Moments calculated with Program I are of varying magnitude but suspiciously 
of only on polarity. This is not in accordance with expectations, namely 
moments varying in both directions. This problem is avoided by cor- 
recting the angles of attack data by a constant value representing the 
difference of the time-averaged values. Several runs of Flight 60 
represent the basis for estimating these correction values. The left 
probe angle of attack is reduced by 2.02 degrees and the center probe 
reduced by 0.67 degrees. Figure 3.3 shows the angles of attack as 
recorded during a 50-second period of Flight 60, Run 18, and Figure 3.4 
shows angles of attack with the corrections included. 

The first direct result of strip theory calculations is shown in 
Figure 3.5. Using the spanwise linear distribution, the lift divided by 
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Figure 3.3. Recorded angles of attack. 
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Figure 3.4. Corrected angles of attack. 







weight obtained in the first calculations with Program I is plotted. 

The lift is normalized with aircraft weight. The aircraft lift^to- 
weight ratio is computed and is compared with a plot of load factor 
recorded by the on-board inertial navigation system (INS) in the figure. 
The theoretical 1 ift-to-weight ratio is computed independently of the 
INS and is based on data collected by use of pitot tube and angle of 
attack probes. Any dynamic influence on lift is neglected; but, since 
the test flight is close to being level and unaccelerated, the dynamic 
terms are expected to be small. The agreement between the two curves 
suggest the strip theory model is reliable in estimating lift on a wing. 
Note that the entire 41.5-second Run 21 of Flight 6 is analyzed, whereas 
only a 50-second period inside the total 140-second Run 18 of Flight 60 
is studied. The first 40 seconds of Run 18 of Flight 60 is omitted as 
is indicated by the time variable of that run in the figures. Also, the 
mean airspeed of Flight 60, Run 18, is 113 m/s but 96 m/s in Flight 6, 
Run 21. With the independent variable as distance (i.e., x = Vt), the 
abscissa would be approximately 15 percent elongated in a plot of 


Flight 60. 

Figure 3.6 represents the l_ ADD in Equation 1.3, which is computed 
as the difference of the two calculations illustrated schematically in 
Figure 3.2. The lift difference is generally- small but can reach values 
up to 10 percent of the weight. The problem with offset due to the 
discrepancy in angle of attack between probes is corrected in Flight 60, 
Run 18, but not in Flight 6, Run 21. This can explain the apparent 
negative time-averaged in Flight 6, Run 21. It is a result of an 
approximately one half of a degree lower time-averaged angle of attack 
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Time (seconds) 


Figure 3.6. Additional lift due to spanwise turbulence 






recorded in the nose of the aircraft, see Figure 3.7. This makes the 
time-averaged lift with uniform distributions greater than the time- 
averaged lift with spanwise distributions which results in the negative 
difference. In this case this difference may simply be caused by upwash 
at the wing tips and not an inclination induced error in the center 
probe. 


Figure 3.8 is identical to Figure 3.6 except that induced drag not 
lift is plotted. is much smaller which is anticipated. 

The roll moment in Equation 1.1 due to spanwise turbulence is 
shown in Figure 3.9. The mean value of in Flight 6, Run 21, is 
negative because of the only slightly higher time-averaged angle of 
attack on the right wing tip compared to the left wing tip. As noted, 
this asymmetry is also displayed in Figure 3.6, and could possibly be 


real because of the relatively short time period over which the calcula- 


tion is made. The roll moment of Flight 60, Run 18, in Figure 3.9 is 
based on the corrected angles of attack as previously mentioned 
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(Figure 3.3, page 23). Similar results for the yaw moment due to span- 
wise turbulence is presented in Figure 3.10. 

The magnitudes of and are of little meaning unless they 
are weighed against known quantities. To provide the reader with a feel 
for the magnitudes of these values, the roll moment developed as a 
function of aileron input and yaw moment due to a rudder deflection are 
plotted in Figure 3.11 . It shows the equivalent to the largest roll 
moment in Figure 3.9 is equivalent to about 10 degrees aileron input. 
Similarly, the yaw moments due to spanwise turbulence can reach magni- 
tudes comparable with a 3 degree rudder deflection. Figure 3.12. 

Figures 3.11 and 3-12 were generated using the aileron control power 

C and the rudder control power C coefficients listed in Appendix B. 

£< 5 a 

The second part of this study investigates the influence on air- 
plane performance generated by the four quantities L A qq, ADD’ ^ADD’ 
and jt add . 
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Figure 3.10. Computed additional yaw moment due to spanwise 
turbulence. 








CHAPTER IV 


FLIGHT SIMULATION USED TO EVALUATE RESPONSE EFFECTS 

In the first part of the study, time histories of aerodynamic 
forces and moments on the wing due to spanwise turbulence are computed. 
In the second part, the effects of spanwise turbulence on the wing are 
incorporated into a 6D0F flight simulation to investigate aircraft 
response. The methods and results of the second part are given in this 
chapter. 


1 . IMPLEMENTING SPANWISE TURBULENCE 

A second computer code is used to implement spanwise turbulence 
into aircraft response analyses. This program is FWG's six-degrees-of- 
freedom flight simulation routine. The wind variation measured during 
a given flight is used as an input to the computer program. The compu- 
ted flight trajectory parameters are then compared to those measured by 
the INS system on the actual airplane when encountering the same wind. 

Two sets of calculations are employed. First a reference case 
which does not include spanwise variations is computed. For the refer- 
ence case only the wind variation measured by the central boom is input 
to the 6D0F program. This wind is assumed to vary with time but to 
act uniformly over the entire aircraft. Next the program is modified 
to include the results of spanwise variation based on the strip theory 
calculations described earlier. In this manner, a comparison of the 
two sets of results can be carried out to assess the influence of span- 
wise variation. Figure 4.1 is a block diagram outline of the proce- 
dure for the aircraft performance analysis. 
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The 6D0F program is the same as used in a study by Frost et al . 
(1984). The original program includes wind gradients in the equations 
of motion (Frost and Bowles 1984), but for this purpose gradients are 
neglected. The only other modification is the incorporation of spanwise 
turbulence. This is achieved by simply adding the four terms (Equations 
1.1 through 1.4) which account for spanwise turbulence into the force and 
moment equations of the 6D0F program. Equations 4.1 through 4.61 
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Only pitching moment in Equation 4.5 remains unaffected by the span- 

wise turbulence, whereas side force in Equation 4.3 is coupled with 

Equation 4.1 and thus influenced by changes in lift due to spanwise 

turbulence. The incorporation of precalculated time histories of span- 

wise turbulence effects into the 6D0F flight simulation is acceptable 

subject to the following condition. The computer simulation must closely 

resemble the real test flight from which spanwise turbulence is calculated. 

In this case, the test flight is a straight level flight with only small 
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pertubations from the trim condition, and thus relatively easy to simu- 
late. Two simulations of the test flight, one including spanwise turbu- 
lence and one without, are performed to reveal the effect of spanwise 
turbulence on aircraft response. 

Note, the following assumption applies to the computer simulation. 
The 6D0F program regards the aircraft as a perfect rigid body. The 
corresponding error introduced is expectedly less for the B-57 than for 
a high aspect ratio wing airplane. The air density and the weight of the 
B-57 is assumed to be constant. Furthermore, the aerodynamics of the 
B-57 are represented by aerodynamic coefficients (Appendix B) that are 
assumed to be linear over the range they are used. The airplane is con- 
trolled through the use of the simple linear gain automatic control 
algorithms. The mathematical control laws are explained in further 
detail in Appendix C. 


2 DATA UTILIZED 

From the B-57 Gust Gradient Program flight data, many test flights 
are available. The immense number of flight parameters recorded during 
a test flight is shown in Table 3.1, page 18 . Early flights, however, 
lacked recordings of pilot control inputs. In the second part of this 
study, such records can be useful for comparison of computed and mea- 
sured results, and consequently Flight 60 was chosen as it became 
available in the spring of 1984. Flight 6 dates back to 1982 and does 
not contain pilot control records and is not analyzed any further. 

In the strip theory calculations, the three angles of attack and 
the corresponding airspeeds are utilized. In the 6D0F simulation program 
the atmospheric conditions measured on Flight 60, Run 18, are used. Air 
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density is calculated from temperature and static pressure. The wind 
components are obtained by removing aircraft motion, measured with the 
INS system, from sensor measurements, the atmospheric inputs including 
wind and density are used to reconstruct conditions of the test flight. 

3. RESULTS OF 6D0F FLIGHT SIMULATIONS 

This section presents the results of the 6D0F flight computer simu- 
lations of Flight 60, Run 18. The two calculations, one incorporating 
spanwise turbulence and one without, are both plotted on the same graph 
to make comparison easy. 

The three wind components of Flight 60, Run 18, as measured with 
the INS during the flight are plotted in Figures 4.2, 4.3, and 4.4, 

Figure 4.3 shows the airplane encountered a strong sidewind a few 
seconds into the time period considered of the test run.* The strong 
changes in north/south wind occurring around the 50-second mark (Figure 
4.4) not surprisingly produce the strongest spanwise turbulence reflected 
in L AD p in Figure 3.6, page 26. 

In Figure 4.5 the recorded yaw angle on Flight 60, Run 18, is 
plotted. About 43 seconds into the run the yaw angle reaches roughly 
negative 11 degrees indicating the nose pointing to the left. This is 
predictable and expected as the airplane entered a strong sidewind from 
the left at that point in time, see Figure 4.3. The recorded airplane 
response is relevant in this context to serve as reference state that 
computer predicted response can be compared with. Refering to Figure 
4.6, it is demonstrated that the computer simulation of Flight 60, 

Run 18, compares favorably with the actual flight. Figure 4.6 also 

*Note that the airplane heading is due north on Flight 60, Run 18. 
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Figure 4.4. Longitudinal wind component (north/south). 


illustrates the effect of adding the spanwise turbulence terms into the 
6D0F flight simulation. With spanwise turbulence included, the amplitude 
of the oscillations generally increase and furthermore closer resemble 
the real flight. 

Figure 4.7 shows the influence of the spanwise turbulence terms on 
the roll response of the 6D0F computer simulation. A noticeable increase 
in roll angle amplitude occurs at the 79- and 85-second marks as a result 
of the spanwise turbulence. Throughout the 50-second simulation, the 
spanwise turbulence excites the higher frequency roll response. 

Figures 4.8 and 4.9 are generated from the 6D0F computer simulation. 
Total roll moment coefficient is given by expressing Equation 4.4 in coeffi- 
cient form. Similarly, total yaw moment coefficient is obtained from the 
coefficient form of Equation 4.6. Again, the difference between the 
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solid and the dashed lines in Figures 4.8 and 4.9 is simply the addition 
of the spanwise turbulence terms of Equations 1.1 through 1.4 into the 
force and moment equations. The increase in high-frequency content of 
Figure 4.8 is due to the term plotted in Figure 3.10, Much of the 
high-frequency is damped out as is evident in Figure 4.7, but the spanwise 
turbulence also contains lower frequencies which affect the roll response 
substantially. In Figure 4.9 the yaw moment coefficient seems to be far 
less affected by the addition of spanwise turbulence terms. This is to 
be expected since Figure 4.6 shows less influence from spanwise turbulence. 
Analyzing the relative influence of the terms in Equation 4.6, it is 
found that the rudder term is responsible for the characteristic shape of 
the yaw coefficient. 

Thus, from the results presented in this chapter, it is demonstrated 
that spanwise turbulence affects the response of the aircraft. The roll 
mode is especially sensitive to variations in turbulence along the wing. 
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Figure 4.6. Calculated yaw angle, 
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Figure 4 9. Calculated total yaw moment coefficient 



CHAPTER V 


CONCLUSIONS 

The results of this study are based on NASA B-57 flight test data. 
The magnitude of spanwise turbulence effects is quantified and also the 
influence of spanwise turbulence on airplane dynamic behavior is 
demonstrated. The severity of the turbulence recorded on the computer 
simulated test flight. Flight 60, Run 18, is best described by the 
2.8 m/s standard deviation of the vertical gust component. 

It is evident from the results that wind variations across the 
B-57 wing can indeed produce large aerodynamic moments on its wings. 
Such roll moments are equivalent of up to 10 degrees of aileron input. 
The influence of spanwise turbulence on the calculated roll mode is 
significant. The pilot workload is also increased. 

The spanwise turbulence also affects the yaw mode of the aircraft 
as is illustrated by the calculations in this study. Since yaw moments 
produced by spanwise turbulence only represent a fraction of similar 
roll moments, the yaw mode is excited, but less severly than the roll 
mode. 

The calculations of aerodynamic forces on the wing showed only 
small changes when spanwise turbulence was replaced by the traditional 
uniform distribution. The greatest difference in lift calculated due 
to spanwise turbulence is 10 percent of the aircraft weight, but gen- 
erally the difference is closer to 2 percent. 

The strip theory lift calculations were expressed as a lift-to- 
weight ratio and compared with the flight data load factor. With the 



gust alleviation factor (see Appendix A) applied to the calculated 
loading, the difference between the two curves is negligible (Figure 
3.5, page 24). The contribution to lift from aircraft dynamics is 
ignored in the theoretical calculations. This comparison suggests 
the strip theory lift calculations are realistic. 

The six-degrees-of-freedom flight simulation is limited by a number 
of necessary assumptions. The aerodynamic modeling of the B-57, in 
particular, is critical. The aerodynamic coefficients and derivatives 
were estimated without any wind tunnel or test flight information. 

With these assumptions made, spanwise turbulence effects are shown to 
substantially influence response. 
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APPENDIX A 


GUST ALLEVIATION FACTOR 

The strip theory calculations of lift in Chapter III are based on an 
instantaneous angle of attack in a steady-state condition. While this is 
approximately true for a gradually applied disturbance, it overestimates 
the loads for a sudden disturbance such as that due to a sharp-edged gust. 
The change of lift is gradual (transient) until steady-state value is 
reached. The response depends upon the characteristics of the gust as 
well as the aircraft. 

For a gradual upgust as illustrated in Figure A.l, the normal accel- 
eration response will be as shown in Figure A. 2 (Babister 1980). The 
time ( t-j ) it takes the airplane to fly through the transition (s-j) is 
related to the airspeed V as follows: 

s 1 = Vt-j (A.l) 

The vertical acceleration increases almost linearly up to time t-| and 
then falls off exponentially. The gust alleviation factor K is defined 
(Babister 1980, p. 108) as the ratio of the maximum vertical accelera- 
tion experienced due to a given gust to the theoretical maximum vertical 
acceleration in a sharp-edged gust with steady-state conditions. As 
expected, the gust alleviation factor decreases as the transition region 
increases. Assuming steady-state conditions, the gust factor becomes 
unity as s^ goes to zero. Furthermore, as the transition region becomes 
large, the unsteady flow effects get smaller. Figure A. 3 which repre- 
sent a plot of the factor K reproduced from Bisplinghoff (1951). Variable 
H is the normalized gust transition region defined as: 
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Vertical Acceleration 




Figure A. 2. Vertical acceleration due to a gradual gust. 
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Steady 



H = s -j / c (A. 2) 

and y is a mass parameter given by: 
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The plot illustrates the importance of the unsteady effects for sharp- 
edged gusts. 

The mass parameter indicates that for an aircraft of low mass, the 
alleviation factor or acceleration ratio gets small as the transition 
region increases. For a high mass aircraft, the acceleration ratio is 
less dependent upon the length of the transition region. However, the 
maximum acceleration itself is inversely proportional to the mass 
parameter and is given by: 
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where m is the maximum vertical gust velocity. 

For the B-57, the following values are obtained under the flight 
conditions in question: 

y g = 24 * 6 
H = 10.6 

K = 0.7 

This value for the gust gradient coefficient is applied to the lift 
calculations in the following way. The difference between lift and 
weight is computed at each time step, and this difference is multiplied 
with K. That fraction is added back to the weight to give the corrected 
lift. The procedure is illustrated in Figure A. 4 for a lift curve which 
is artificially constructed. 
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Figure A. 4. Correction to lift using the gust alleviation factor K. 


APPENDIX B 


B-57 CHARACTERISTICS USED IN SIMULATION OF FLIGHT 60, RUN 18 
1. TRIM CONDITIONS FOR FLIGHT 60, RUN 18 


Trim airspeed 
Air density 
Tail incidence angle 
Thrust force 
Aircraft mass 

Level flight, no flaps, gear up 

2. AIRPLANE 

Wing area 
Wing span 

Mean aerodynamic chord 
Aspect ratio 
Moments of inertia 


v = 113.0 m/s (370.7 ft/sec) 
p = 0.96186 kg/m 3 (0.00187 slugs/ft 3 ) 
i' t = -0.22 degrees 
T = 36,509 Newtons (8,207 Ibf) 

W = 20,420 kg (45,017 lb) 

PARAMETERS 

S = 89.2 m 2 (960 ft 2 ) 

b = 19.5 m (63.95 ft) 

c = 4.57 m (15 ft) 
mac 

AR = 4.3 

I = 325,300 kg - m 2 (240,000 slug ft 2 ) 

XX 

I = 271,100 kg - m 2 (200,000 slug ft 2 ) 

I = 650,600 kg - m 2 (480,000 slug ft 2 ) 

I = 11,000 kg - m 2 (8,100 slug ft 2 ) 

L. A 
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3. AERODYNAMIC COEFFICIENTS AND STABILITY DERIVATIVES 

C. = 0.07 
0 

C L = 4.98/rad 

a 

C, = 0.419/rad 
6 e 

C, = 0.670/ rad 
1 t 

C, - 2.90/rad 

q 

C^ # = 1.45/rad 

a 

C D = C. 2 /13.5 
o L 

C m * 0.064 
m o 

C m = -1.0/ rad 

m a 

C m = -0.860/ rad 
6 e 

C m = -1.72/rad 

\ 

C m = -9.65/rad 

m q 

C m = -5. 83/ rad 

m* 

a 

C = -0.286/rad 

y g 

C = -0.071 /rad 
y x 
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c = 0.057 3/ rad 

n e 

c = -0.0321 /rad 

n 6 

r 

C = -0.00286/rad 

C n = -0.070/rad 

C = -0.200/rad 
n r 

C = -0.0286/rad 
C = 0.0040/rad 

V 

C. = 0.0458/rad 

\ 

C = -0.5/ rad 
P 

C = 0.081 /rad 

X/ - 
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APPENDIX C 


CONTROL LAWS 

The s i x-degrees-of -freedom simulation uses the following automatic 
control laws: 


= DEI (V, - V Q ) + DE2(e - e trim ) 


<$ a = DAl(V y - V g ) + DA2 • <j> 


V = o R Hv y - V f 


ref 

ref 

ref 


) + DR2 • 


where: 

DEI = -0.0025 (Elevator control gain factor #1 ) 

V e = 2.0 m/s 

ref 

DE2 =0.06 (Elevator control gain factor #2) 

9 trim = 1,203 de 9 rees 
and 

DAI = -0.005 (Aileron control gain factor #1) 

V a = -0.08 V w 

a ref y 

DA2 = -0.80 (Aileron control gain factor #2) 
and 

DR1 = 0.015 (Rudder control gain factor #1) 


= -0.12 V. 


(C.l) 

(C.2) 

(C.3) 


ref 


DR2 =0.50 (Rudder control gain factor #2) 
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APPENDIX D 


FLIGHT INFORMATION 

1. Flight 6, Run 21 

Date: July 14, 1982 

Location: Denver, Colorado 

Project: Joint Weather Airport Studies (JAWS) 

Start time: 14:41:58 MDT 

Duration: 41.5 seconds 

Altitude above ground: 1000 feet 

Mean airspeed: 96 m/s (315.0 ft/ sec) 

Flight description: Level flight 

2. Flight 60, Run 18 

Date: February 1, 1984 

Location: Boulder, Colorado 

Project: Orographic Effects Campaign 

Start time: 16:57:15 MST 

Duration: 140 seconds 

Altitude above ground: 1500 feet 

Mean airspeed: 113 m/s (370.7 ft/sec) 

Flight description: Level flight 
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